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High energy vibrational scattering in the binary systems UC and US is measured using time-of-flight inelastic 
neutron scattering. A clear set of well-defined peaks equally separated in energy is observed in UC, correspond¬ 
ing to harmonic oscillations of the light C atoms in a cage of heavy U atoms. The scattering is much weaker in 
US and only a few oscillator peaks are visible. We show how the difference between the materials can be under¬ 
stood by considering the neutron scattering lengths and masses of the lighter atoms. Monte Carlo ray tracing is 
used to simulate the scattering, with near quantitative agreement with the data in UC, and some differences with 
US. The possibility of observing anharmonicity and anisotropy in the potentials of the light atoms is investigated 
in UC. Overall the observed data is well accounted for by considering each light atom as a single atom isotropic 
quantum harmonic oscillator. 

PACS numbers: 63.20.dd, 63.20.Pw, 78.70.Nx 


I. INTRODUCTION 

The complex electronic structure of uranium leads to a wide 
variety of unusual and diverse behavior in the uranium salts 
UX (X = C, N, P, S, As, Se, Sb, Te, Bi), and therefore the 
materials in this family have been the subject of many past 
investigations characterizing their magnetic and vibrational 
properties ifTHTOl . The phonon spectra in the UX systems 
have been comprehensively studied, however there has been 
renewed interest in more detailed calculations of their vibra¬ 
tional properties as certain members of the family, including 
UN and UC, are under active consideration for next genera¬ 
tion nuclear fuels ifTTUT^ . 

Recently, new and unexpected features were discovered in 
the vibrational spectrum of UN ITtII . The same single crys¬ 
tal of UN studied in earlier work||7l [15] [TbI was re-examined 
via modem neutron time-of-flight (TOE) spectroscopy using 
the SEQU0IA|[T7] [TSl and ARCS llT9ll spectrometers at the 
Spallation Neutron Source, Oak Ridge National Laboratory. 
Speciflcally strong vibrational scattering was found at ener¬ 
gies above the usual acoustic and optic phonon branches. A 
series of evenly-spaced, high energy modes was observed, and 
detailed quantitative analysis showed that these modes could 
be attributed to the nitrogen atoms behaving as independent, 
isotropic, 3D quantum harmonic oscillators (OHOsl lfTHI^ . 
Some of the features of the QHO modes, for example the in¬ 
trinsic broadening, were consistent with predictions of a bi¬ 
nary solid model (211. This model applies to systems with two 
different types of atoms with disparate masses and explains 
how the motion of the heavy atoms affects the QHO behavior 
of the lighter ones. 

The observation of these well-deflned modes at energies 


above the highest optic phonon in an ordered single crystal 
contrasts with the conventional view of vibrational response 
in crystalline solids. Above the highest optic phonon modes 
the response is usually weak and relatively featureless 1221 . 
Prior to the UN discovery, exceptions to this rule were gen¬ 
erally found in binary metallic hydrides ll23H27l . for example 
in ZrHa; systems, where hydrogen atoms occupy interstitial 
sites 1^ . However in those systems, the hydrogen modes 
usually exhibit signiflcant anisotropic and anharmonic effects, 
mainly due to H-H interactions, crystalline anisotropy, and 
the diffusion of the H atoms ll23ti26l . Typically only a few 
modes are observed. Conversely, the nitrogen oscillations in 
UN show well-deflned peaks up to the 10^^ order. 

Known pre-requisite conditions for the QHO modes to be 
clearly observable in binary alloys include a large mass ratio 
between the light and heavy atoms, and weak interactions be¬ 
tween light atoms. Although the aforementioned binary solid 
model provides some guidance, there is little experimental in¬ 
formation on the dependence of this part of the vibrational re¬ 
sponse function on factors such as the mass ratio and atomic 
neutron cross-sections. In this work, we investigate how these 
factors affect the high energy vibrational scattering in UX via 
time-of-flight neutron spectroscopy measurements on single 
crystals of UC and US. Eor UC, we And a series of well- 
deflned high energy vibrational modes analogous to those ob¬ 
served for UN. Overall these modes are described well by 
a QHO model for the carbon atoms. Conversely, while the 
US data also shows evidence for the high energy vibrational 
modes, only the lowest few are observed and they are much 
weaker in intensity compared to those seen in UC and UN. 
We discuss the reasons for these differences. In addition we 
explore the effect of multiple scattering on the observed spec¬ 
trum, and the possibility of extracting information on the di- 
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rectional dependence and anharmonicity of the light atom po¬ 
tential via measurements of the QHO modes. 


II. INELASTIC NEUTRON SCATTERING 

All neutron scattering measurements reported here were 
collected using the SEQUOIAO \M and ARCS (191 TOP 
Fermi chopper spectrometers at the Spallation Neutron Source 
of Oak Ridge National Laboratory. The same depleted ura¬ 
nium single crystals of UC and US used in previous studies 
were investigated in this work ll28ll^ . Both samples had sim¬ 
ilar total volumes on the order of 1 cm^. For the neutron scat¬ 
tering experiments, each single crystal was mounted in an alu¬ 
minum can and loaded in a closed cycle helium refrigerator. 
All data were collected at T = 4 K, with the [HHL] scattering 
plane horizontal. A Fermi chopper was used to obtain several 
different incident neutron energies, including Ei = 80 meV, 
250 meV, 500 meV, 700 meV and 800 meV. The details for 
each chopper setting are given in Table |I| All datasets were 
normalized against a vanadium standard to account for vari¬ 
ations of the detector response and the solid angle coverage. 
Empty can measurements were performed at T = 4 K and 
subtracted from the US datasets, however a small amount 
of A1 that made up the sample mount was not perfectly ac¬ 
counted for. We discuss below the implications of this imper¬ 
fect background subtraction on the scattering observed in US. 
A similar background subtraction was found to be unneces¬ 
sary for the UC data because the intrinsic features of interest 
are intense. 



[-2 -2 L] 


FIG. 1: The phonons along the [-2 -2 L] direction for (a) UC (ARCS) 
and (b) US (SEQUOIA) with Ei = 80 meV. In both cases, the acous¬ 
tic and optic phonon modes are well-separated. The optic modes’ 
signal to background ratio is much higher for UC as compared to 
US. The energies of the optic modes for UC and US are known to be 
48 meV and 40 meV Q. The white text O and A indicate the general 
location for the optic and acoustic phonon modes, while B denotes 
significant background contribution. 


TABLE I: List of experimental conditions used to collect data for (a): 
UC, and (b): US, on SEQUOIAiL;^ = 80 meV only) and ARCS(Ei = 
80, 250, 500, 700 and 800 meV). Here dpc is the slit spacing of the 
Fermi chopper, Rfc is the radius of curvature of the Fermi chopper, 
ufc is the frequency of the Fermi chopper, uto is the frequency 
of the TO chopper, and AE FWHM represents the corresponding 
instrumental energy resolution at zero energy transfer. 


(a) UC, T ^ 4 K 


Ei (meV) 

80 

500 

700 

800 

dFC (mm) 

1.5 

0.5 

0.5 

0.5 

Rfc (m) 

0.58 

1.53 

1.53 

1.53 

ufc (Hz) 

120 

480 

600 

600 

uto (Hz) 

90 

180 

180 

180 

AE FWHM (meV) 

9.0 

18.2 

26.1 

31.4 


(b) US, T ^ 4 K 


Ei (meV) 

80 

250 

500 

dFC (mm) 

3.6 

0.5 

0.5 

Rfc (m) 

1.53 

1.53 

1.53 

ufc (Hz) 

240 

360 

480 

UTO (Hz) 

90 

180 

180 

AE FWHM (meV) 

4.1 

8.8 

18.2 


Figure shows representative phonon measurements with 
Ei = 80 meV for (a) UC from ARCS and (b) US from SE¬ 
QUOIA along the [221/] direction in reciprocal space. The 


results for the acoustic and optic phonons are consistent with 
those previously reported in the literature (71 EH. Both the 
acoustic and optic modes in UC are clearly visible with large 
signal compared to background. On the other hand, the optic 
phonon modes of US in particular are only weakly visible in 
the data. The decreased intensity of these excitations can be 
understood by considering the one phonon structure factor, g, 
for the rocksalt structure. Along symmetry directions this can 
be written as (TOllSOl : 

= (1) 


with 


/ bydu ^ bxdx \ 

V \/M a/to ) 


( 2 ) 


where Q is the neutron momentum transfer, ^ is a unit vector 
describing the phonon polarization, M is the mass of U, m is 
the mass of the light atom X, du and dx are the relative dis¬ 
placements of the U and the X ions, subject to the constraint 
E d?x — and bx are the neutron bound coher¬ 

ent scattering lengths (TOl . The sign between the two terms in 
Eq. (2) depends on the (HKL) indices of the Brillouin zone, 
where (+) and (-) correspond to all even and all odd indices 
respectively. 
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For optical phonons, at the zone center \duM\ = \dxTn\ 
and therefore g' can be rewritten as: 


_ rnM ^ bu ^ bx ^ 


(3) 


In the limit of M 


oo : 


/)2 

hm ^ 

M^oo m 


(4) 


This expression indicates that the intensity of the optic 
modes is determined primarily by the mass and scattering 
length of the light atoms. Under the assumption of similar 
scattering lengths, the intensity of the optic modes decreases 
with increasing mass m of the light atom. Table [n| further em¬ 
phasizes this point by comparing several parameters relevant 
to the phonon intensities of UC, UN, and US EH. It is seen 
that the optic phonon cross section for US is more than an or¬ 
der of magnitude smaller than the comparable cross-sections 
for UC or UN. A1 is also included in the table since. It leads 
to significant background scattering in the experiments, as dis¬ 
cussed below. 


TABLE II: Relevant parameters for phonon intensities in UC, UN, 
US, and Al. (H] 


Atom 

m (amu) 

b(fm) 

b^/m {fm^ / amu) 

N 

14 

9.4 

6.31 

C 

12 

6.6 

3.63 

s 

32 

2.8 

0.25 

u 

238 

8.4 

- 

Al 

27 

3.4 

- 


was incorporated into the fit. The decaying ex¬ 
ponential term is expected from a simple diffusive model of 
multiphonon scattering |[32j|33l. are the mode positions 
that were fit independently and the last term incorporates all 
others sources of background. 

The solid line in figure represents the fitted curve from 
Eq. (5). The fitted peak positions of the modes for all mea¬ 
sured datasets are shown in Table and are evenly spaced. 
The average spacing between the modes, incorporating all 
fitted datasets, was found to be hujQ = 48±1 meV for UC. 
The values are averaged by a weighted approach meaning 
hujQ = (ficci -|- E Tiujji/Ti) jTi. Values 

calculated from individual peaks, i.e.ficcn/^ deviate from the 
average by less than 1 meV, which implies a highly harmonic 
potential. 


TABLE III: Peak positions of the high energy vibrational modes (in 
meV) for UC from fits to the Q-integrated data. Please note that the 
error bars are statistical from the fitting of the data and do not account 
for instrumental resolution or systematic effects. 


(a) UC, T ^ 4 K 


n\Ei 

500meV 

700meV 

SOOmeV 

1 

47.7(3) 

- 

- 

2 

96.2(3) 

98.1(2) 

99.7(4) 

3 

142.8(4) 

144.1(3) 

146.0(5) 

4 

189.1(6) 

190.9(3) 

193.0(4) 

5 

235.6(1) 

237.3(4) 

239.0(4) 

6 

283.8(36) 

283.7(6) 

284.7(6) 

7 

- 

329.9(7) 

331.0(1.0) 

8 

- 

372.9(1.2) 

375.5(1.7) 


average Uujq. 48±1 meV 


Figure depicts the orientationally averaged, high energy 
response observed for both UC and US on ARCS. Several 
evenly-spaced vibrational modes are easily visible in the UC 
Ei = 500 meV spectrum shown in Fig.j^a). In the US data 
with Ei = 500 meV (Fig. [^b)) the high energy vibrational 
modes are less obvious. On the other hand, some weak modes 
can be ascertained in the higher resolution = 250 meV data 
depicted in Fig. [^c). In both materials the modes appear on 
inspection to be evenly spaced. 

We first discuss the scattering seen in UC. Figureshows 
the Ei = 500 meV Q-integrated data for UC, with two dif¬ 
ferent panels corresponding to linear (top) and log (bottom) 
y-axes. Modes up to 7th order are clearly visible as peaks in 
the data. 

To test whether these are evenly spaced the Q-integrated 
dataset is fitted to the following functional form: 


I{E) = Y,Tne 


{E-EnV 

9/-r2 


+ Be 


^-XbE 


+ ^0 


(5) 


Tn is a scale factor for each individual Gaussian peak, and 
is the standard deviation for each peak. An empirical back¬ 
ground consisting of a constant Bq and decaying exponential 


The fact that the spectrum is characterized by a series of 
evenly-spaced modes suggests that the light atoms behave like 
independent 3D QHOs similar to UN lfTHl . This hypothesis 
can be tested quantitatively by comparing the relative intensi¬ 
ties of the modes to the known dynamical structure factor for 
the mode of a QHO at low T, which has a simple analyti¬ 
cal form given by: 


Bfl 


D f fiQ^ 

n\ \2muJo 


n 

exp 


(zEE) 
\ 2muJo J 




( 6 ) 


where E{uj) = S{huj — nhujo) and D is a constant. 

Assuming that one has a carbon atom QHO with the mea¬ 
sured energy spacing of hujo = 48 meV and a mass of 12.01 
a.m.u. the calculated value of hl2mujQ is 0.0036 A^. Given 
the range of Q measured and the Q dependence in equa¬ 
tion (6), it is surprising that the higher order modes are vis¬ 
ible in figure 3. This can be understood by considering that 
Sn{Q^oj) only accounts for single scattering events. With 
equally spaced modes, however, the contributions from mul¬ 
tiple scattering will also peak at energies corresponding to 
mode positions, and as discussed below the multiple scatter¬ 
ing contributes to intensity at low Q. 

We consider the meaning of a single scattering event in 
both QHO language and an alternate, but equivalent descrip- 
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FIG. 2: Color contour plots of the experimental neutron scattering spectra for (a) UC, Ei — 500 meV, (b) US, Ei — 500 meV and (c) US, 
Ei = 250 meV. The evenly-spaced, high energy vibrational modes are clearly visible in the UC plot, but much less pronounced in the case of 
US. The small Q inelastic scattering seen in US arises from magnetic excitations. 



FIG. 3: The Q-integrated intensity versus energy plot for UC at 4 K, 
with Ei — 500 meV. The upper plot shows the systematic decline of 
intensity for higher modes. The lower (log) plot shows that the QHO 
modes remain visible for this dataset up to the 7th mode. The black 
points are from the experimental data and the solid red line is the fit 
to the data using Eq. (5). 


tion based on creating Einstein phonons of fixed frequency. 
At T = 0 an inelastic neutron scattering event, at an energy of 
nfiujQ in QHO language, represents a transition from the QHO 
ground state to the eigenstate. In Einstein mode language 
it consists of a single scattering event, with total wavevector 
transfer magnitude Q, that creates n Einstein phonons, each 
of energy hujo. In either case the scattering cross-section con¬ 
tains a factor On the other hand, a multiple scattering 


event in QHO language corresponds to a set of multiple tran¬ 
sitions from the ground state to the excited state, while 
in the phonon description several different scatterings occur 
creating in aggregate n Einstein phonons. The scatter¬ 
ing event contributes a factor Q‘j to the total observed cross- 
section. These processes are depicted in Eigurej^for the case 
n = 2. The single scattering event shown in Eigure|^(a) will 
have an intensity I ^ and a requirement that the vector 
sum of the individual Einstein phonon momenta Qi and Q 2 
has a magnitude equal to Q. On the other hand, in the multiple 
scattering event shown in Eigure (b) each created Einstein 
phonon contributes a separate factor to the measured inten¬ 
sity which is therefore proportional to {QiQ 2 )^. The total net 
vector momentum transfer Q measured in the experiment can 
take on a wide range of values, including Q = 0 as depicted 
in the figure. The average over all possible multiple scattering 
combinations can lead to a cross-section that is independent of 
Q. Therefore the observed intensity of the higher order peaks 
at small Q arises principally from multiple scattering. 

With this in mind, we fit the constant-cuts of the UC data 
to a modified expression of 5'(Q, cc) for QHOs ifTdll . given by: 

SniQ, w) = A^Q‘^^exp{-CQ‘^) + (7) 

Here C is ?il2muj{) and a Q-independent Bn term is included 
to account for the multiple scattering described above. Eor 
the ideal QHO model, AnjD = jul. By relaxing this con¬ 
straint and allowing An and C to be independent parameters 
in the data fitting, one can gain a sense of how much the UC 
data deviates from the ideal QHO limit. 

Eigure |^a) depicts constant-cuts from the 500 meV 
dataset for the n = 1-5 modes, centered about nficco for the 
nth (20 meV integration range). Similar cuts were made 
for the Ei = 700 meV and 800 meV datasets (not shown). 
Eigurej^b) plots the ratios of the fitted parameters AnjAi for 
the Ei = 500 meV, 700 meV, and 800 meV datasets. The 
solid line indicates the prediction for the QHO model with 
Tiujq = 48 meV, with the non-integer n values interpolated 
by using r(n + 1) to calculate n\. As in the case of UN, 
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(a) single scattering event 




(b) nnulti-scattering event 
•32 



FIG. 4: Schematics showing the difference between (a) single scat¬ 
tering process and (b) multiple scattering processes leading to the 
observed neutron scattering spectra for UC, UN, and US. ki and kf 
refer to the incident and final neutron momenta. In (b) values of Qi 
and Q 2 were chosen such that the total momentum transfer Q = 0 
but the scattering intensity is non-zero (see text for details). 


there is excellent agreement between the QHO model and 
the data over 20 orders of magnitude. The experimental C 
value, corresponding to the zero point motion of the oscil¬ 
lator, also agrees well with the model. By incorporating all 
three datasets into the fit, the average value of C was found 
to be 0.0031(1) A^, which is close to the calculated value for 
the ideal QHO of 0.0036 discussed above. These findings 
provide strong confirmation that the high energy vibrational 
modes observed in UC correspond to quantum oscillations of 
the carbon atoms in the system. 

Figure shows the Q-integrated data for US with 
Ei = 25{) meV and Ei = 500 meV, including fitted curves 
using Eq. (5). The same fitting approach, described earlier for 
UC, was also used for the US data. However due to the weaker 
mode intensities in US, it is not possible to fit as many peaks 
in the US dataset. Nonetheless the modes are evenly spaced, 
similar to UC, with an averaged hujQ of 41(1) meV. The fitted 


peak position of each individual mode is shown in Table IV 


TABLE IV: Peak positions of the high energy vibrational modes (in 
meV) for US from fits to the Q-integrated data. Please note that the 
error bars listed are from the fitting of the data and do not account 
for instrumental resolution. 


US, T ^ 4 K 


n\Ei 

250meV 

500meV 

1 

40.9(3) 

- 

2 

79.9(9) 

85.6(2.4) 

3 

- 

122.9(3.6) 

4 

- 

166.7(5.8) 


average hujo : 41 zb 1 meV 


The weak intensities of higher QHO modes in US can be 
largely explained by the phonon structure factor of the rock- 
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FIG. 5: (a) The Q-dependence of the intensity for the n = 1-5 

UC oscillator modes with Ei = 500 meV, integrated over the energy 
ranges shown in the legend. The solid lines are fits to the QHO model 
described in the text, while the (+) symbols represent constant- 
cuts from the MC simulations presented in section III. (b) The An 
coefficients for the n = 1-5 modes normalized by the values of Ai. 
The coefficients are extracted from fits to the Ei = 500 meV (black 
triangles), 700 meV (red circles) and 800 meV (blue squares) ARCS 
datasets. The solid line is the prediction of the QHO model with 
huo = 48 meV and m corresponding to the mass of a carbon atom. 
The ratio is plotted on a logarithmic scale and spans almost 20 orders 
of magnitude. 


salt structure, as discussed earlier by Eq. (3) and (4), and Ta- 
ble|n| In short the greater mass and shorter scattering length of 
sulfur results in a much smaller b‘^/m compared to carbon or 
nitrogen. At small Q the scattering at the nth 

mode position is 

dominated by multiple scattering processes sequentially creat¬ 
ing single optic phonons. The observed cross-section is there¬ 
fore proportional to Eor this reason the higher QHO 

modes appear much weaker in US compared to UC. 

The correction for the Al background scattering is also a 
significant complication for the observation of QHO modes in 
US. The QHO cross-section for each mode reaches a maxi¬ 
mum at the value of Q where the mode energy corresponds to 
the energy of the recoil scattering, which is given for an atom 
of mass m by ErecoU = Ti^Q^/2m. As seen in Table II the 
masses of Al and S are relatively close, so the recoil scatter¬ 
ing from the Al is strong near the position where one expects 
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E (m eV)' 

FIG. 6: The Q-integrated intensity versus energy plot for US at 4 
K, with Ei = 250 meV and Ei = 500 meV. The QHO modes are 
much less intense compared to those observed in UC, with only 2 
modes observable for Ei = 250 meV and 4 for Ei = 500 meV. 
The black points are from the experimental data and the solid red 
lines are the fits to the data using Eq. (5). Please note that the Q- 
integrated range for the US data is Q > 5 only. This is done to 
avoid contributions from the magnetic excitations at low Q as seen 
in Fig. 2 (b) and (c). 

to see the scattering due to the QHO modes of S. Altogether, 
quantitative analysis of the intensity of the QHO modes in US 
is more difficult. 


III. MONTE CARLO RAY TRACING SIMULATIONS 

Monte Carlo (MC) ray tracing simulations of both the 
UC and US ARCS neutron scattering experiments were per¬ 
formed to gain a better understanding of the various fac¬ 
tors contributing to the observed scattering intensity. This 
approach worked exceptionally well to describe the neu¬ 
tron scattering spectra observed for UN l^ . The ability to 
model more realistic instrument and sample configurations 
is ever increasing 1341 [SSj. The simulations used the MC 
ray tracing program Monte Carlo Virtual Neutron Experiment 
(MCViNE), which was based off Mcstas|[36l [37l and devel¬ 
oped in the Distributed Data Analysis for Neutron Scattering 
Experiments (DANSE) (381 software development project. In 
each simulation, the sample was modeled as a cube with a 
1 cm^ volume, approximating the shape of the sample used 
in the measurement. The configuration files for the simula¬ 
tions were created through a series of simple modifications to 


the files used for the UN MC simulations described in detail 
in Ref. (20l, allowing for a straightforward extension of those 
calculations to UC and US. 

Intrinsic broadening of the QHO modes is included in the 
simulation, as it is observed experimentally and predicted by 
the binary solid model (211 . The broadening arises from the 
fact that the heavy U atoms are not completely stationary, and 
can be modeled by replacing the Dirac-Delta function in the 
QHO expression with a Gaussian: 

N f (huj — nhuJo)‘^\ 

- ) ( 8 ) 

where the Gaussian width r(T) is a function of temperature 
and given by: 

M here is the mass of the heavy atom, in this case the U atom. 
Z{u) is the acoustic phonon density of states calculated with 
a Born-Van Karmen model and 0 is the maximum band fre¬ 
quency for the acoustic phonon modes. 

A Monte Carlo ray tracing simulation for UC with 
Ei = 500 meV is shown as a color contour plot in Eig. [TJa) 
and the Q-integrated result is presented in Eig. [^d). This 
simulation is analogous to the one reported in Ref. (^ for 
UN and includes various sample kernels accounting for elas¬ 
tic scattering, QHO scattering (from carbon atoms), acoustic 
phonon scattering, and all forms of multiple scattering arising 
from the processes described above. Eigure [^a) described 
earlier includes cuts from the Ei = 500 meV MC simulation 
as (+) symbols with the experimental data, where C is fixed at 
the value of 0.0036 corresponding to a carbon QHO with 
Tiujq = 48 meV and a mass of 12.01 a.m.u. 

Monte Carlo simulations for US, incorporating sample ker¬ 
nels with the same ingredients as for UC, are shown as color 
contour plots in Eig. Bb) and (c) with intrinsic broaden¬ 
ing included. The Q-integrated results are also presented in 
Eig. [Tj^e) and (f). Eor (f), data simulated without intrinsic 
broadening is also included. At low signal levels the mea¬ 
sured data approaches a constant value, above the simulations. 
This indicates that there is a background contribution, limit¬ 
ing the measurement sensitivity, that is not captured in the 
Monte Carlo simulations, nor by the A1 can subtraction. The 
additional background could arise from a combination of ef¬ 
fects including imperfect A1 background subtraction, multiple 
scattering involving the magnetic response or events partially 
external to the sample, and possibly a contribution from fis¬ 
sion neutrons in the sample that moderate in the instrument 
shielding. 

When comparing the high energy vibrational modes of UC, 
UN, and US it is apparent that the differences in the U:X mass 
ratio (U:C = 19.8, U:N = 17, U:S = 7.4) lead to more obvious 
QHO behavior in the lighter atoms. As the U:X mass ratio 
decreases, the dispersion of the optic phonon modes has been 
shown to increase systematically (T0]| . shifting the yibrational 
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FIG. 7: Color contour plots from the Monte Carlo ray tracing simulations described in text are shown here for (a) UC, Ei = 500 meV, (b) US, 
Ei = 500 meV and (c) US, Ei = 250 meV. The simulations that incorporate intrinsic broadening of the QHO modes agree well with the data 
for UC, but not for the case of US. (d), (e), and (f) show the corresponding Q-integrated plots for the three experimental data sets, and compare 
them directly to the results from the simulations, with solid green lines showing the fitted curve for each of the experimental data using Eq. (5). 
(f) also shows MC simulation results both including and excluding intrinsic broadening of the QHO modes. 


behavior of the light atoms further away from the localized 
QHO picture. A smaller U:X mass ratio has also the effect of 
pushing the optic phonon frequencies down closer in energy 
to the acoustic modes, so any resulting QHO excitations tend 
to be more closely spaced and therefore harder to resolve for 
a given instrumental energy resolution. 


IV. ANHARMONICITY AND ANISOTROPY 

All QHO data presented to this point for both UC and 
UN lfm have been powder-averaged, so only Q is considered 
and not Q. The reason this approximation works so well is 
because the modes are highly isotropic. However, in princi¬ 
ple the spacings of the modes could reflect both anisotropy 
and anharmonicity in the local potential. To see whether these 
effects are experimentally observable in the QHO modes of 
UC, we collected Bl Ei = 700 meV dataset with very high 
statistics along three crystallographic directions, [100], [110], 
and [111]. The C atoms have nearest neighbor distances 
along these directions in the ratio of l:v^:v/3, as one can 
see from considering the crystal structure shown in Fig.j^a). 
The effects of anisotropy would be manifested by different 


mode spacings for the different directions, while anharmonic¬ 
ity would be apparent by different mode spacings as a func¬ 
tion of energy. Unfortunately, both of these features will be 
masked by the effect of multiple scattering, since at any value 
of Q it tends to cause peaks to appear at locations correspond¬ 
ing to nhujQ. As a result, a straightforward plotting of oscilla¬ 
tor peaks along different directions does not show any signifi¬ 
cant anisotropy or anharmonic effects, and additional analysis 
is needed to account for the multiple scattering. 

As an effort to minimize the effect of multiple scattering 
present in the single crystal data, we take advantage of the 
fact discussed earlier that the small Q scattering for the higher 
modes is totally dominated by multiple scattering, and use the 
low Q data as a pseudo-background, as described here: For 
each of the three directions we integrate over a range perpen¬ 
dicular to that direction given by = [-2,2] A“^. As a 
signal we use the high Q part of the scattering defined by 
20 < Qii <30 A“^, and as a background we utilize the re¬ 
gion 0 < Qii <8 A“^. The subtraction of these two leads to 
constant-Q difference cuts along the three directions as shown 
in Fig. [^b). Note that the n = 1 peak is not sufficiently re¬ 
solved in the Ei =700 meV data to be included here. 

The difference cuts are fit to Eq. (5) to determine mode po- 
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FIG. 8: (a) The crystal structure of the binary uranium systems UX. The light X atom is in an octahedral cage of U atoms. The colored lines 
indicate the three primary crystallographic directions, [001], [110] and [111]. The carbon atom nearest neighbor distances in these directions 
are in the ratio of l:\/2:\/3. (b) Constant-difference cuts along the major crystallographic directions [001], [110], and [111] for UC. The 
solid curves are fits of the data to Eq. (5). (c) Peak positions along the three different directions taken from the fits in (b). The dotted line 
corresponds to the expected peak positions for an isotropic QHO. (d) The difference between the peak positions of the modes from our UC 
data and a perfectly harmonic potential, (e) Local potential of the carbon atoms calculated by DFT. The solid curves are fits to the anharmonic 
function V(x) given in the text, (f) The difference between the peak positions of the UC data and the result from the DFT calculations. 


TABLE V: Peak positions of the QHO modes (in meV) extracted 
from fitting difference cuts of the Ei = 700 meV single crystal UC 
data along three major crystallographic directions. The form of the 
difference cuts is explained in the text. 


n 

[001] 

[110] 

[111] 

1 

2 

102(3) 

98(3) 

99(3) 

3 

147(2) 

142(3) 

144(3) 

4 

192(2) 

192(2) 

191(2) 

5 

239(3) 

239(2) 

235(2) 

6 

296(3) 

287(2) 

282(3) 


sitions along the [001], [110], and [111] directions; fits are in¬ 
dicated in Fig.j^b) by the solid curves, and the fitted peak po¬ 
sitions are listed in Table|V] These are also shown in Fig.j^c), 
with the mode positions determined from the data and 

the dashed curve representing the expectation for an ideal car¬ 
bon QHO with Tiujq = 48 ±1 meV. Fig.j^d) shows the differ¬ 
ence between the mode positions for the ideal model (nhujo) 
and the experiment. The error bars are statistical and are larger 
than those for fits to the powder averaged data, presumably 
due to the lower number of counts in this single crystal dataset. 


From the fits, the average oscillator energies are given for the 
different directions as F^ooi=49.0(6)meV, F^iio=48.0(6)meV, 
and F^iii=47.9(6)meV. There is no compelling evidence for 
anisotropy in the potential although the experimental data 
along the [001] direction did exhibit a slightly higher fre¬ 
quency than those along the other directions. This observation 
is compatible with density functional theory (DFT) calcula¬ 
tions discussed below. 

To gain insight into any expected anisotropy and anhar- 
monicity of the oscillator modes in UC, we have used DFT to 
calculate the potentials seen by the C atoms (see the appendix 
for details). The calculated potentials for the three principal 
directions are shown by the symbols in Fig.j^e). The poten¬ 
tial in each direction is described well by the following ex¬ 
pression: 

V{x) = ax^ -h px^ (10) 

with the results shown by the solid lines in the figure. The val¬ 
ues determined for a and /3 are also indicated. The calculation 
shows a small amount of anisotropy, and also indicates that the 
anharmonicity should be maximized along the [001] direction. 
The latter result suggests that the harmonic approximation be¬ 
comes less valid for in the direction with the shortest nearest 






















9 


neighbor distance. Indeed the fitted P for [001] is more than 
twice the magnitude of the /3s along the other directions.. 

Assuming that the term is small, the energy levels corre¬ 
sponding to this potential can be calculated according to first 
order perturbation theory (391, with the energy level cor¬ 
responding to: 

/ ^ \ 2 

Er^ = hujQ^n + 1/2) + ( -- ) + fin + 3) (11) 

\lmujQ J 

and the energy difference between the level and the n = 0 
ground state given by: 

En — Eq = nhujQ + 6^C‘^{n‘^ + n) (12) 

with C = hl2muJo. 

Assuming that the first mode occurs at 48 meV, the ener¬ 
gies of higher QHO levels along each of the 3 directions can 
be calculated using Eq. (12) along with each direction’s cor¬ 
responding a and /3 shown in Fig.j^e). Fig.j^f) displays the 
difference between the oscillator mode positions determined 
from the experiment and the value E^^^ from first order per¬ 
turbation theory. The downward slope indicates that the ex¬ 
perimentally determined anharmonicity is not as large as the 
estimate. One must bear in mind that despite the attempt to 
correct for multiple scattering in the data, it is greater at higher 
order peaks and therefore could still mask what is a rather sub¬ 
tle effect. 

V. CONCLUSIONS 

For the binary crystal UC, time-of-fiight neutron scattering 
measurements reveal a series of well-defined, equally spaced, 
high energy vibrational modes that can be attributed to quan¬ 
tum harmonic oscillator behavior of the carbon atoms in this 
system. Measurements of the QHO modes along the high- 
symmetry [001], [110], and [111] directions reveal that these 
excitations are characterized by at most only a very small 
amount of anharmonicity and anisotropy. Similar time-of- 
fiight neutron scattering data for US also shows evidence for 
the high energy vibrational modes, but only a few are clearly 
observed and they are much weaker in intensity. The differ¬ 
ence can be understood by considering the U:X mass ratio, 
the quantity jm that characterizes the scattering strength of 
the modes, and the interference of A1 and S recoil scattering 
with the signal in US. Some progress has been made in mod¬ 
eling the various contributions to the scattering, but a better 
way to handle multiple scattering will be needed to accurately 
measure anharmonicity and anisotropy in the QHO modes. 
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APPENDIX: DET CALCULATION 

For UC, the electronic structure within the density func¬ 
tional theory (DET) was obtained using the Quantum 
ESPRESSO package HO). The calculation was performed 
using a plane-wave basis set and an ultrasoft pseudopoten¬ 
tial (411 optimized in a RRKJ scheme (421. The uranium pseu¬ 
dopotential was obtained from an ionized electronic configu¬ 
ration: with cutoff radii equal to 3.5 atomic 

units (a.u.), 1.7 a.u., 2.6 a.u. and 1.6 a.u. for d and / 
angular momentum. The electronic levels deviate from the 
all-electron ones by less than 0.1 meV. We used the Perdew, 
Burke, Ernzerhof (^ exchange-correlation functional. The 
Brillouin zone (BZ) summations were carried out over a 
4x4x4 supercell. The electronic smearing with a width 
of 0.02 Ry was applied according to the Methfessel-Paxton 
method. The plane wave energy and charge density cut-offs 
were 73 Ry and 1054 Ry respectively, corresponding to a cal¬ 
culation accuracy of 0.2 mRy/atom. The carbon atom po¬ 
tential was obtained from the total energy modification of a 
2x2x2 supercell when one carbon atom was shifted from 
the equilibrium position in the [100], [110], or [111] direc¬ 
tions and the remaining atoms were held fixed in their equi¬ 
librium positions. 
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